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We have studied the plasmonic properties of aperiodic arrays of identical nanoparticles (NPs) 
formed by two opposite and equal graded-chains (a chain where interactions change gradually). We 
found that these arrays concentrate the external electromagnetic fields even in the long wavelength 
limit. The phenomenon was understood by identifying the system with an effective cavity where 
plasmonics excitations are trapped between effective band edges, resulting from the change of pass- 
band with NP’s position. Dependence of excitation concentration on several system’s parameter was 
also assessed. This includes, different gradings as well as NP’s couplings, damping, and resonant 
frequencies. In the spirit of the scaling laws in condensed matter physics, we developed a theory 
that allows us to rationalize all these system’s parameters into universal curves. The theory is quite 
general and can also be used on many other situations (different arrays for example). Additionally, 
we also provided an analytical solution, in the tight-binding limit, for the plasmonic response of 
homogeneous linear chains of NPs illuminated by a plane wave. Our results can find applications 
on sensing, near field imaging, plasmon-enhanced photodetectors, as well as to increase solar cell 
efficiency. 

PACS numbers: 73.20.Mf, 73.21.-b, 78.67.-n, 42.79.Ek 


INTRODUCTION. 

Concentration or focusing of electromagnetic fields 
(EMFs) has been a goal of great interest almost since 
the inventions of lenses. This is due to its potential ap¬ 
plications, which in modern times span from imaging and 
sensing to lasing. Traditionally mirrors, lenses, or com¬ 
binations of them, such as resonant cavities 111,1 , have 
been used for this purpose. But the diffractive nature of 
electromagnetic waves has always imposed a hard limit to 
these devices. There have been several proposals to over¬ 
came this, including super lenses made of negative refrac¬ 
tive index metamaterials JU Q, superoscillations [5], and 
time reversal focusing [7]. However, difficulties associ¬ 
ated with them, such as complex designs involving active 
materials or the need for an almost absolute control of 
the incident EMF, have undermined their development 
into concrete applications. 

Probably the most successful strategy for subwave¬ 
length, or even deep-subwavelenght, concentration of 
light is the exploitation of evanescent waves. All though 
predicted a long time ago [1, 0], we had to wait until the 
end of the 20th century, when we acquired enough control 
of matter at nanoscale, to see those predictions turn into 
experimental results J0-14]. Despite its success, there 
are however certain drawbacks in using evanescent waves 
for light concentration, as for example the volume-versus- 
confinement problem. 

It is known that certain nanostructure presents under 
illumination regions of high EMFs or hotspots. Now, 
if one want to increase even further the EMFs inside 
a hotspot one can sharpen the tips or the region of 
the nanoparticle (NP) where concentration takes places. 


However, this also reduces the volume of the hotspots 
0, an effect that can not be neglected if the target 
molecule to be sensed is a macromolecule for example, or 
the goal is to increase the efficiency of solar cells [16| . 

Another possibility, involves the use of arrays of NPs 
or nanoscopic structures where radiation concentration 
occurs inside the system as consequence of the interac¬ 
tion among its different components 0. Here, two dif¬ 
ferent situations should be distinguished. The first one 
occurs when nanostructures have one or more dimen¬ 
sions approaching the excitation wavelength and thus 
it is necessary to consider retardation effects. The con¬ 
cepts behind many of these retardation-based plasmonic 
light concentrators are based on scaled radio frequency 
antenna designs 01 [l7l-[20j. There are of course other 
alternatives, such as the use of aperiodic metallic waveg¬ 
uide arrays [2l|. The second situation occurs when the 


size of the nanostructure is significantly smaller than the 
wavelength of the incident light. Here, the entire struc¬ 
ture experiences a uniform electric field at any instant 
and then the quasistatic approximation is valid. Exam¬ 
ples of this are given by the use of fractal aggregates 
[22[ or self-similar chains HU of subwavelength particles 
usually called nanolenses. It is worth noting that in the 
mentioned and in many other examples, the lack of peri¬ 
odicity of the arrays seems to be key in EMF’s focusing. 

Aperiodic arrays of NPs, sometimes called graded ar¬ 
rays if the parameters describing the array vary in a 
graded way, have shown interesting and uniques features. 
Such as the appearance of new types of excitations called 
gradons [24], or a continuous frequency-depend localiza¬ 
tion of excitations in different part of the system HU- 

In this work we study a kind of aperiodic array of NPs, 
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FIG. 1. Scheme of the system studied, a mirrored graded- 
chain of NPs. It consists of two opposite and equal graded- 
chain of identical metallic NPs. 


that we called mirrored graded-chain of NPs. They are 
formed by two equal but opposite graded-chains of NPS, 
see Fig. [l] Array’s geometry allows it to concentrate 
external EMFs, but with the important property that 
excitation concentration does not come at the expense of 
reducing hotspot’s volume, as in the case of nanolenses. 
Through the work we focus on the physics behind the 
phenomenon and the effect of different system’s param¬ 
eters on excitation concentration. In the spirit of the 
scaling laws in condensed matter physics, we develop the 
theory that allows one to build universal curves for the 
properties studied. The theory, valid also for other ar¬ 
rays of NPs, requires only arrays made of equal NPs that 
can be treated in the quasistatic limit. 

The work is organized as follows: In section I we briefly 
introduce the method used for the calculations. In sec¬ 
tion II we develop various aspects of the behavior of lin¬ 
ear chains of equidistant NPs that will result useful in the 
next sections. In section III we present the main results 
of our work. This section is subdivided into two main 
subsections “Nearest-neighbor approximation” and “De¬ 
pendence on system’s parameters “. In the first one we 
use a simple version of our system to extract the physics 
behind excitation concentration. While in the second 
one, we systematically study different aspects that could 
affect that. In this last subsection we develop a simple 
theory that allows us to build universal curves for the 
frequency of the maximum excitation concentration and 
the strength of it. There is also a short subsection at 
the end ’’Comparison with simple graded-chains” where 
we briefly compare excitation concentration between mir¬ 
rored and simple graded-chains. Finally, in section IV we 
summarize the main conclusions. 
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where £?t° tal = E^ xt) + & n,m, eo is the free space 

permittivity, V is the volume, e is the dielectric constant 
of NP’s material, e me d is the dielectric constant of the 
host medium, and L is a geometric factor that depends 
on NP’s shape and the direction of E relative to the array 
0. 

For a linear array of NPs, transversal and longitudi¬ 
nal excitations, T and L respectively, do not mix, which 
allows us to write F n?m as [3ll. l32l|: 
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where d n?m is the distance between NPs, and the complex 
constant 7 depends on the wavenumber k of the excita¬ 
tion and the orientation of the NP’s array relative to the 
direction of P. 


7 n,m = I 1 - ikd n,m ~ {kd n , m ) 2 }e lkdn ’ m 

ln,m = -2[! - ikd ntm \e ikdn ’ m . (3) 

Note that in the quasistatic approximation, kd « 0, 
and 7 l are 1 and —2 respectively. 

The dipolar moments P n , and the external electric 
fields E^ xt can be arranged as vectors P and E resulting 
in: 

P = (Iw 2 - M) _1 RE = xE, (4) 


where x is the response function, M is the dynamical 
matrix and R is the diagonal matrix that rescales the 
external applied field according to the shape, volume and 
material of the NP. 

Assuming a particular model for e in Eq. [T| one can 
obtain the expressions for the different elements of M and 
R. For ellipsoidal NPs and using a Drude-Sommerfeld’s 
like model for e, e = eoo — c<; 2 /(u ; 2 + iujrf), one obtains 
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COUPLED DIPOLE APPROXIMATION. 
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In this work we study linear arrays of NPs where the 
interactions among them are modeled using the coupled 
dipole approximation 26|-32]. In this model, each n th - 
NP is described by a dipole moment P n induced by the 
electric field produced by the other dipoles m at n, P n ,m, 
and the external electric field at n, E ^ 1 . Assuming a 
generic ellipsoidal shape for the NPs, this gives for the 


and 


Rn,n — e of^n^p( n )/- (7) 

where cj p ( n ) is the plasmon frequency of the NP with in¬ 
dex “n”, ^ sp ( n ) is its resonant frequency, r n (cj) is its 
decay rate, and uj^ n n is the coupling between NPs 
n and m. The factor / can in principle depend on a;, 
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FIG. 2. Typical dispersion relation for longitudinal and 
transversal excitations of NP’s homogeneous-linear-chains. 
Calculations were done assuming only nearest-neighbor in¬ 
teractions, r) = 0, and a quasistatic approximation. 


for Coo 7 ^ e me d, which typically sli ghtl y shifts the reso¬ 
nances of the system, see Ref. |3l|, [32|. However, as we 
are mainly interested in extracting general trends of our 
system, we will take / as constant for simplicity. For 
the same reason, we will consider only the case of equal 
NPs with the same resonant frequency for transversal 
and longitudinal excitations, and a linear form for r(cj), 
i.e. T (cj) = rjuj, which neglects radiation damping. 

According to the above expressions, the role of cj sp is in 
general to set the scale of the problem, which is reflected 
as a shift in the spectra or in the dispersion relations. 
For this reason and in order to gain generality, from now 
on we will only use rescaled variables for the excitation 
frequencies, coupling constants, damping terms, Ru , and 
the self energy n (see next section) of the linear arrays. 
That is, those variables should be understood as c J — 
u/u SP , SP > v' = v/u SP , R'u = Rii/u 2 P , and 

n' = n/u;^ respectively. 


HOMOGENEOUS LINEAR CHAIN OF 
NANOPARTICLES. 


In this section we discuss different aspects of the behav¬ 
ior of homogeneous linear chains of NPs, a system where 
distances among neighboring NPs are equal. The discus¬ 
sions here presented will result useful in the next sections. 
In order to obtain simple analytical expressions that will 
allow us to understand our problem in simple terms, we 
will treat our systems in the quasistatic approximation 
and neglect interactions beyond first neighbors. 

For an homogeneous linear chain, the dipolar mo¬ 
ment of two NPs can only differ in a phase factor, thus 
Pn = P m e _ ^ m_n ^ d where d is the distance between 
neighboring NPs. Considering only nearest-neighbor in¬ 
teractions and negligible damping terms, one can find the 
value of the wave vector q for a given frequency of exci¬ 


tation, i.e. the dispersion relation of our system 00, 

u 2 (q ) = 1 + 2 uj 2 cos (qd). (8) 

Fig. [2] shows an example of this, where it can be seen 
that changing the type of plasmonic excitation not only 
shrinks the passband, the interval of frequencies for which 
Eq. [8] has solution, but it also inverts the curve’s shape 
due to a change in the sign of c The maximum and 
minimum values of the passband are known as band edges 
and for this particularly simple example are given by: 

w|(<?) = l±2w 2 . (9) 

Now, let us analyze the excitation of a linear array of 
NPs illuminated by a plane wave. The electric field of 
the external source at NP “n” is 

E^ xt) = E 0 Ee ink * d (10) 

where Eq is the modulus of P, E its direction, n = 
(..., —1, 0,1, 2...) the NP’s index, and k z d = kdcos(0) 
with 6 being the angle between the NP’s chain and the 
incident light. In a homogeneous linear chain of NPs 
with only nearest-neighbor interaction one can calculate 
analytically all the elements of the response function x, 
defined in Eq. [4| jib, [H[ 
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where all the self energies n^’s are equal in the limit of 
infinite NPs ( N -T oc) 
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Therefore, when the whole chain is illuminated, the dipo¬ 
lar moment P n of a given NP can be calculated from: 
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where Eq ' Ms the transversal or one of the longitudinal 
components of EqE. 


Fig. [3] shows an example of a linear homogeneous array 
of NPs illuminated perpendicular to the chain. As can 
be seen, the analytical solution overlaps the numerical 
results. The maximum of each spectrum occurs at one 
of the band edges, the lower one for longitudinal excita¬ 
tions, negative cj 2 , and the higher one for transversal ex¬ 
citations, positive uj 2 . This behavior can be understood 
in terms of Eq. [IT] in the limit of small 77 . For example 


in the case of longitudinal excitations the fraction 



is exactly “ 1 ” at the lower band edge and exactly “- 1 ” at 
the higher band edge. That means that two consecutive 
terms of Eq. [14] contribute constructively in the former 
case and destructively in the latter one (the other way 
around for transversal excitations). In more physiccal 
terms, that means that the absorption of light by two 
consecutive nanoparticles interfere contructively or de¬ 
structively if the excitation frequency is one or the other 
band edge. Furthermore, the density of states has peaks 
only at the band edges, so only there it is expected a con¬ 
siderable absorption of the external field. In Fig. [3j we 
can also notice that the height of the peaks for T and L 
excitations are approximately the same as that of an iso¬ 
lated NP, note that the normalization factor is precisely 
the maximum value of |P | 2 of an isolated NP. The small 
differences in the peaks’ high are due to the frequency de¬ 
pendence of the damping terms. This is also reasonable 
as all NPs are being excited and absorb in the same way 
and the resulting excitation must be distributed without 
a preferential direction through the array. 


Finally, it is useful to briefly discuss the effect of sev¬ 
eral deviations from the ideal model considered above. 
As can be found in several references (26l432| . the effect 
of taking into account all orders of interparticle couplings 
is essentially a widening of the passband plus a shifting of 
its center. Retardation effects on the other hand, have a 
stronger influence on the properties of the chain, driving 
the appearance of new resonances and valleys product 
of constructive and destructive interferences in the inter¬ 
particle couplings. Those effects are stronger at higher 
frequencies and when the size of the problem (NP’s radii 
or interparticle distances) is not negligible compared with 
the wavelength of the excitation. However, in this work 
we will focus on systems with the appropriate scale such 
as to safely neglect these retardation effects. At the end 
of the last section we will retake this issue. The presence 
of defects in the chain, such as NPs of different radii 
or inhomogeneous NP’s separations, induces the appear¬ 
ance of localized states, which correspond to collective 
excitations of a reduced number a NPs whose resonance 
generally lies outside the passband [ 3 Q|. 


(M 



FIG. 3. Square dipolar moment, |P| 2 , of a NP in a 
homogeneous-linear-chain as function of the square frequency 
of the external field, uj 2 . The wave vector of the external 
field is perpendicular to the chain with E longitudinal (L) or 
transversal (T) to it. The nearest-neighbor approximation as 
well as the quasistatic approximation are considered, 77 = 0.01 
and uj 2 = —0.1 and 0.05 for L and T excitations respectively. 
All curves are normalized to the maximum value of |P | 2 of an 
isolated NP with 77 = 0.01. Numerical results correspond to 
|P | 2 of the central NP of a finite array of 81 NPs. 


MIRRORED GRADED-CHAIN OF 
NANOPARTICLES. 


Now that we have discussed some relevant aspects of 
linear arrays of equidistant NPs, we can focus on under¬ 
standing the effect of changing from homogeneous inter¬ 
actions to slightly graded couplings. As shown in Fig. [TJ 
the system studied consists of an inhomogeneous linear 
array of metallic NPs where the distance between neigh¬ 
boring NPs incresases from the center to the ends of the 
array while the size and shape of NPs are kept constant. 
Two different models for distance’s increasing are studied 

d n = doa\ n \ model 1 (Exponential) , , 

d n = do(1 + \n\/3) model 2 (Linear) 

where a and /3 are the incremental factors, n = 0 cor¬ 
responds to the central NP, and we have simplified the 
notation for d by replacing d n?n+ 1 (or d_ n? _( n+1 )) by d n 
(or d- n ). 

We will see that this type of systems focus the external 
electromagnetic field on the central NPs. In what follows 
we will first analize our system in the nearest-neighbor 
approximation. Then, we will disscus the effect of taking 
into account different aspects that approach our simple 
model to a more realistic situation, but always trying to 
find general trends in its behavior that can result useful 
for designing potential applications. 
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FIG. 4. \P\ 2 vs uj 2 for different NPs in a mirrored graded- 
chain of NPs (see inset). Nps’ distances are changed according 
to model 1 (Eq. fl6l) with a = 1.035. The system is illumi¬ 
nated perpendicularly to the array with E longitudinal to the 
chain. Nearest-neighbor approximation as well as the qua¬ 
sistatic approximation are considered. We used uj 2 — —0.1, 
77 = 0.01, and the total number of NPs is 81. All curves are 
normalized to the maximum value of \P\ 2 of an isolated NP 
with 77 — 0 . 01 . 


Nearest-neighbor approximation. 


Fig. 0] shows an example of the “spectrum” (where 
spectrum stands for \P\ 2 vs uj 2 ) of the studied system. 
As we can see, now the maximum of the “spectrum “ of 
the central NP is slightly shifted from the band edge po¬ 
sition toward the center of the passband. New resonances 
appear, especially for NPs far from the center, but more 
interestingly, there is a considerable increase of the peak’s 
height of the central NPs together with a decrease of it for 
NPs away from the center. Indeed, not only the height 
of the spectrum is increased but also its total area. 

Fig. [5]shows a 3D color map of \P\ 2 vs uj 2 and NP’s po¬ 
sition. The plot presents a typical seagull’s contour with 
well resolved resonances around the middle that slowly 
turns into a continuum toward the edges. The general 
shape of the plot can be understood by following the ap¬ 
proach of treating the system as “locally” homogeneous 
|2f[ l24j | . That means that even though the properties of 
the array change continuously with NP’s position, these 
changes are smooth enough to treat the system locally as 
an homogeneous linear array. Thus, we can define a po¬ 
sition dependent band edge at every NP’s position. The 
continuous green lines of Fig. [5] mark exactly that. Note, 
the almost perfect matching between the green line and 
the maximum of the excitation for NPs away from the 
center. 


Xiao et.al. [24] studied infinite linear arrays of NPs 
and harmonic oscillators with graded interactions in a 
given direction, which would be equivalent to take only 
one half of our system (see Fig. [1]). The key difference 


FIG. 5. (Color online) - 3D plot of \P\ 2 vs uj 2 and NP’s 
position in the array. Conditions are the same as in Fig. [4j 
Continuous green lines show the band edges of the local equiv¬ 
alent homogeneous chains (see text and Eq. Grey dashed 
lines mark the three lowest eigenvalues of M for longitudinal 
excitations. 


between their and our case is that here an excitation at a 
given frequency can not travel indefinitely in a given di¬ 
rection, even with 77 = 0. For example, take an excitation 
that start at the crossing between the grey dashed line at 
the highest frequency and the descending green line that 
correspond to the lowest passband edge (the lower left 
one). This excitation can travel through the chain only 
to the right (propagation to the left will put the excita¬ 
tion outside the passband). This excitation can continue 
traveling but only until it reaches the other crossing with 
the band edge, the one at the right. Beyond that point 
the excitation would be again outside the passband. This 
implies that excitations are necessarily trapped between 
the crossing with the band edges and thus the system 
acts as an effective optical cavity. 

As it is well known, optical cavities can increase enor¬ 
mously electromagnetic fields inside them due to coherent 
accumulation of excitations. This explains not only the 
concentration of excitations in our system but also the 
discrete nature of resonances around the center of the 
system. Obviously, due to finite damping, excitations 
that has to travel distances greater than a characteristic 
decay length to interfere with themselves will not present 
this “optical cavity” effect. This explains why the system 
behaves as homogeneous linear chains for NPs far apart 
from the center. 

Calculation of cavity’s resonant frequencies is not that 
simple in our system as the effective cavity changes its 
length with frequency and its group velocity with posi¬ 
tion. However one can always calculate numerically the 
eigenvalues of M. Grey lines in Fig. [5] show the real part 
of the three lowest eigenvalues of M. Those are the val¬ 
ues at which the resonant condition occurs. Naturally 
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n (index of the NP) 


FIG. 6. The three first normal modes of a mirrored graded- 
chain of NPs. Conditions are the same as in Fig. [4j 



not all of them are bright modes as the excitation field, 
a plane wave perpendicular to the array, can not excite 
anti-symmetric modes for example. The eigenmodes cor¬ 
responding to those eigenvalues are shown in Fig. [6l Note 
that the lowest eigenmode has, as expected, its maximum 
at the center of the array which explains why the high¬ 
est concentration of radiation occurs around the central 
NP. For transversal excitation, modes are similar but P n 
changes its sign between neighboring sites. The spectra 
are also similar as well as the 3D plots but, as expected, 
the peaks for the different NPs appear close to their high¬ 
est band edges (with |^(L)/^I = I Co ’x(t)D- Besides that, 
the only difference is that peak’s heights are lower due to 
a higher damping (T = rju). 




FIG. 7. Upper figure - The same as Fig. [5]but a — 1.01 and 
all interactions are considered. Lower left figure - Vertical 
cuts of the upper figure at n = 0, 3, and 40 (same as Fig. 
[4] but under the present conditions). Lower right figure - 
Horizontal cut of the upper figure at uj 2 = 0.78 (response of 
the system at constant frequency of excitation). 


Dependence on system’s parameters. 

In this section, we will consider different factors that 
deviate the system from the simple model discussed pre¬ 
viously. In all the figures shown, we have taken |Po|^ ax , 
the maximum value of |P| 2 of the central NP, as an in¬ 
dicator of excitation concentration. Another possibility 
could have been, for example, to take the area under the 
curve of the spectra of |Po| 2 vs uj 2 . However, as the re¬ 
sults proved to be qualitatively the same, we will not 
discuss the last one. 

Fig. 0 shows an example of the response of a graded- 
chain of NPs under the same condition as in Fig. [5] but 
considering all interactions among NPs, not only first 
neighbors. Qualitatively the behavior of the system is 
essentially the same. The differences are that the concen¬ 
tration of radiation is slightly increased and the peak’s 
positions are shifted towards lower frequencies, both are 
better appreciated in the lower left sub-figure. The shift¬ 
ing is expected as the band edges of NP’s waveguides are 
red-shifted when considering all interactions among NPs 
and longitudinal excitations [261432] . Apart from that, the 
only difference is that NPs of the end of the chain (NP’s 
index -40 and 40 in the plot) do not show yet the response 
of isolated NPs as in Fig. [5] This is also reasonable as 


the effective interaction between these NPs and the rest 
of chain is larger when we consider interaction beyond 
first neighbors. Thus, a larger array would be required 
to observe that. 

The lower right sub-figure of Fig. 0 shows the square 
dipolar moment of different NPs of the chain for a fixed 
frequency of excitation. Note, how quickly excitation 
concentration decays when moving away from the central 
NPs. This effect is interesting because the system, which 
is a 1-D planar structure, behaves similarly to the tips of 
a SNOM though 0- A phenomenon that, in principle, 
can be used for near field spectroscopies. 

Fig. [8] is an example of the convergence of |Po|^ ax , 
the maximum value of |P| 2 of the central NP, with the 
total number of NPs in the array. Note that a graded- 
chain of only 20 or 40 NPs, depending on cj 2 , shows a 
better concentration of external excitation than any finite 
or infinite chain of equidistant NPs. Interestingly, both 
models for graded-chains, see Eq. [T6l yield almost the 
same result for the values of a and /3 used. Those values 
correspond roughly to the optimal ones (see Fig. [9j). The 
reason for that is simple. Even though a graded-chain of 
NPs is always better than a chain of equidistant NPs, 
the optimal one corresponds surprisingly to an array of 
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FIG. 8 . Maximum value of |P | 2 of the central NP of arrays 
with different length (TV). Three values of co ’ 2 are used: -0.25, 
-0.1 and -0.05 with 77 = 0.01. blue continuous lines, mirrored 
graded-chains where distances are incremented exponentially 
with a equal to 1.005, 1.02, and 1.02. Red dotted lines, the 
same but distances are incremented linearly with /3 = 1 — a. 
Black dashed lines, homogeneous linear chains of NPs. Lowest 
curves correspond to smallest a; 2 . 

almost equidistant NPs, i.e. a and /3 close to 1 and 0 
respectively. Expanding d n /do around a = 1 for the 
exponential model, gives 

a n « 1 + n(a — 1) + 0(a 2 ). (17) 

Thus, 

dn (p) | m odell ~ = OL - l)| model2 (18) 

For that reason, hereafter we will exclusively discuss the 
results for d n given by one of them, the exponential 
model. 

Let us now discuss the effect of a on excitation con¬ 
centration. It is clear from Fig. [8] that there should be 
an optimal value of a between 1 and 00, as a = 1 corre¬ 
spond to equidistant NPs and a = 00 to an isolated NP. 
Fig. [9] shows precisely that. However, the optimal value 
of a is very close to 1 as discussed. Indeed, the bigger the 
coupling (u; 2 ), the closer to 1 results this optimal value. 

The position of the peak of |Po(cj)| 2 , cj^ ax , also changes 
with a. At a = 1, it corresponds to the position of 
the lower band edge of the equivalent homogeneous NP’s 
chain. At a = oc it is 1, which corresponds to the value 
of an isolated NP. Fig. flOl shows this behavior. As can be 
seen in the figure different couplings give different curves. 
However, all curves can be rescaled to a single universal 
one by means of a simple linear transformation, see Fig. 
EH This behavior can be understood by analyzing how 
the eigenvalues of the M matrix, Eq. 01 change with a. 
First, we have to rewrite the matrix M as: 


(19) 


FIG. 9. Maximum value of |P | 2 of the central NP of a mir¬ 
rored graded-chain as function of a (distances change expo¬ 
nentially with “n”). Calculations are converged with respect 
to N and 77 = 0.01. 


where the matrix M', defined through its elements as 
n = 0 and m = (do/d n;m ) 3 , depends only on the 
geometry of the array. For example for d n given by the 

exponential model, = 1/ (^i^ 1 . Clearly a 

matrix U that diagonalize M also diagonalize M' which 
implies: 

ED = uj 2 W + (1 — ir]u:)I (20) 

where U -1 MU = ED and EJ _1 M / U = ED', being ED and ID' 
the diagonal matrices of the eigenvalues of M and M' 
respectively. Therefore, the eigenvalues of M which give 
the position of system’s response maxima are given by 
the eigenvalues of a matrix that only depends on the 
geometry of the array (a in our case) and cj 2 through 

a ) = ^x^o/tvOt) + 1 ( 21 ) 

where KDq^ is the lowest eigenvalue of M' for longitudinal 
excitations or the highest one for transversal excitations, 
and it depends only on a in our case. 

Similarly, one can find the universal behavior of |Po| 2 
by noticing that the eigenvectors of M also depend only 
on the geometry of the array. Using Eq. 0J we can write 

Pn = E Un ^ ~2 - ,) U^Rm,m E{ * Xt) (22) 

, u- - Ui,i ’ 


Now, as all NPs are the same P m , m = Ro and the system 
can be treated in the quasistatic approximation = 

then 

(23) 

uj 2 - (oj 2 D' u (a) + 1 - m ? w ) 

By using this equation one can in principle calculate any 
property of interest of our system by knowing its geomet¬ 
ric contributions, U and (£>'), and the rescale parameters 


Pn _ 
RoEo , 


M = cj 2 M / + (l-ir(cj))I 
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FIG. 10. Position of the spectrum’s peak of |P | 2 vs co 2 , ^ ax) 
for the central NP of a mirrored graded-chain as function of a. 
Calculations are converged with respect to N and 77 = 0.01. 


FIG. 12 . Same as Fig. [9]but using a rescaled maximum square 
dipolar moment of the central NP. 



FIG. 11. Same as Fig. [TOlbut using a rescaled square resonant 
frequency of the central NP. 


Po, an d V which do depend on the material, size and 
shape of individual NPs. Under certain conditions the 
expression can be simplified even further. For example, 
one can assume that the main contribution to the above 
summation at uo 2 = co^ ax is given by l = 0/N. Then the 
following relation holds 


|P 0 (^x,7?,a)|max 

\Ro\ 2 \E 0 \ 2 


Uo,0/N U 0/N, 


w max 


(24) 


where the numerator on the right hand side depends only 
on the geometry of the array (a in our case). Fig. [12] 
shows the validity of this equation for different values of 
cj x and a. We can see that the convergence is excellent 
especially for large values of a or high couplings. 

By interpreting the system as an effective cavity, we 
can easily understand the deviations of Eq. [24] For (a)s 
close to 1 , the curves that represent the band edges as 
function of NP’s positions, see for example the continuous 
green line of Fig. [5j becomes flater which will make the 


eigenvalues of M to collapse. The curvature of these lines 
depend on both cj x and a. The smaller the a and/or cj x 
the less pronounced the curvatures and the less separated 
the eigenvalues of M. Therefore, it is reasonable to see 
deviation of Eq. [24] for a close to 1 and small cj x . 

Fig. [T3l shows an example of the variation of |Po|^ ax 
with 77 which gives an almost perfect l/rj 2 dependence as 
expected. At large ( 77 ) 3 , compared with the separation 
among eigenfrequencies, the approximation in Eq. [24] 
obviously breaks down which is the origin of the small 
deviations at large ( 77 )s. 

It is interesting to highlight that Eqs. [24] and [ 21 ] (or 
some equivalent) are valid beyond the particular model 
treated here of mirrored graded-chains. In particular Eq. 
I2T1 (or the equivalent one for other P^s and its corre¬ 
sponding Dijs) is valid for any system of equal NPs in 
the quasistatic approximation, while Eq. [24] also requires 
values of 77 smaller than the eigenfrequencies’ separation, 
i.e. well defined resonances. Retardation effects affects B 
and the numerator of Eq. [24] Thus, in a sense, they mod¬ 
ify the geometry or the effective geometry of the array. 
Regretfully this turns the clean linear relation between 
Bq^ and cj^ax i n t° a complex non-linear problem as 
now B' = B'(a, fcodo?^), see Eq. [26] Anyway, provided 
the necessary conditions are met, Eqs. l2llandl24l(or even 
Eq. [23] in the case of not well resolved resonances) can be 
used to find universalities in the behavior of any system 
of interest even if we can not solve it analytically. 

Returning to the particular system studied here, it is 
important to address the question of how sensitive is our 
proposed structure to defects in the array. Fig. fT~4l shows 
how | Po | 2 changes when imperfections are introduced into 
the array. In the upper figure we changed randomly the 
distances between NPs (d^ ew = d° ld + Ad) by using a nor¬ 
mal distribution to sample Ad/d° ld . In the lower figure 
we did the same but changing NP’s resonance frequency. 
The width of the normal distribution used to sample the 
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FIG. 13. Maximum square dipolar moment of the central NP 
of a mirrored graded-chain as function of 77 . a = 1.02 and 
cj 2 = —0.1. Green line is a fit to a 77“ 2 function. 

distances (or resonant frequencies) is what is indicated 
in the x axis as “rmsd of d n ” (“rmsd of cjgp, n ”). We 
run 10 calculations for every condition and calculate the 
average value of |Po| 2 and its standard deviation, which 
is indicated as an error bar in the plots. We can see that 
the system exhibits a good fault tolerance. For up to 3% 
or 5% of noise the energy concentration remains almost 
the same. Being, the fault tolerance to defects in the 
distances (d n ) better than that to the resonant frequen¬ 
cies (cjgp(n)). This tolerance to defects approaches our 
proposal to its implementation. Furthermore, it is im¬ 
portant to highlight that when energy concentration is 
based on retardation effects, systems result, in general, 
very sensitive to defects in its components. 

Finally, in order to study retardation effects but with¬ 
out resorting to particular examples, with a given value 
of d, go>sp, etc..., we define fco as. 

&o=^sp/c (25) 

where c is the velocity of light, and rewrite kd in Eq. [3] 
as 

kd = uj(d n:rn /do)(kodo) (26) 

where lj is as usual the renormalized excitation frequency 
(the true uo divided by cj sp ) and kodo is the parameter 
that will account for retardation effects. Fig. [15] shows 
an example of the consequences of increasing retardation 
effects. If kodo is still small, we see a plateau and then a 
smooth decreasing of energy concentration on the central 
NP. When retardation effects increase even further there 
are, as expected, very narrow resonances at some precise 
values of kodo that produce a high concentration of the 
external excitations. However, the resonances are so sen¬ 
sitive to the specific value of fco and do that any defect 
in the distances or a shift in the frequency of excitation 
would destroy it. Furthermore, at a fixed frequency of 



rmsd of d n (%) 



rmsd of (%) 

FIG. 14. Upper figure -Average value and root-mean-square 
deviation (rmsd) of |Po|max f° r mirrored graded-chains with 
random defects in d n . We used a normal probability distribu¬ 
tion to sample the random defects in d n . The X-axis shows 
the rmsd of d n and the height of the vertical bars the rmsd 
of the values of |Po | 2 obtained, a = 1 . 02 , uj 2 = — 0 . 1 ,and 
77 = 0.01. - Lower figure - The same but considering ran¬ 
dom defects on the resonant frequency of individual NPs. 

excitation, increasing kodo is equivalent to enlarge the 
whole system which is contrary to our goal, to develop 
nanostructures that focus external fields in the smallest 
possible region. 

Comparison with simple graded-chains. 

Finally, in this subsection we compared the excitation 
concentration of mirrored graded-chains with that of sim¬ 
ple graded chains. A simple graded-chain corresponds to 
one half of the system discussed above, equivalent to con¬ 
sider only NPs with indexes from 0 to (TV — l)/2 in Fig. 
[1] for example. Fig. [TGl shows excitation concentration in 
terms of the maximum value of \Pi \ 2 vs NP’s positions 
for a sufficiently large chain. The figure can be compared 
with Fig. [8] for large TVs as the same three pairs of oj 2 
and a values where used. As expected, simple graded- 
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k 0 do 


FIG. 15. Maximum square dipolar moment as function of re¬ 
tardation effects (see Eq. EH) for the central NP of a mirrored 
graded-chain with a — 1.005, = —0.25,77 = 0.01. 



n (index of the NP) 


FIG. 16. Maximum square dipolar moment as function of 
NP’s position for simple graded-chains. The total number of 
NPs is 41 and 77 = 0.01. 

chains also behave as resonant cavities. However, they 
present lower excitation concentrations relative to that 
of mirrored graded-chains, especially for high values. 
Furthermore, the position of excitation’s maximum along 
the chains depends on the interaction strength, which 
makes them harder to control. 

CONCLUSIONS. 

We have studied a kind of nanostructure not previously 
reported, up to our knowledge, that we called mirrored 
graded-chain of NPs. We showed that they are able to 
concentrate or focus light even when system’s sizes are 
negligible compared with the wavelength of the external 
excitation. The phenomenon was understood by inter¬ 
preting the system as an effective cavity where plasmon- 
ics excitations are trapped between effective band edges, 
consequence of the change of the passband with NP’s 
positions. 


Dependence of excitation concentration on several sys¬ 
tem parameters was also assessed. Including, different 
models for distances changing, the degree of changing of 
distances with position, as well as NP’s couplings, damp¬ 
ing, and resonant frequencies. A method to build uni¬ 
versal curves for the property of interest was developed 
to understand all these system’s parameters. The key 
idea is based on splitting system’s response into geomet¬ 
ric components, that in general need to be evaluated nu¬ 
merically, and other system’s parameters which will be 
used to rescale the property of interest. The method is 
quite general and can be used on many other situations 
as in principle it only requires arrays made of identical 
particles (nano- or not) well described in the quasistatic 
limit. Note also the similarity between Eq. [4] and the re¬ 
sulting equations for a set of classical coupled harmonic 
oscillators j[34| or the dynamic of a single particle in quan¬ 
tum mechanics [Ip, EliJ . Anyway, provided the mentioned 
conditions are fulfilled, the proposed method can result 
especially useful for situations where system’s complexity 
limits the possibility of analytical solutions but general 
trends were anyhow necessary. 

There are a number of applications where mirrored 
graded-chains of NPs can be useful. In the case of sen¬ 
sors, it can help to solve the problem of hotspot’s volume 
vs intensity. It is known that electromagnetic fields be¬ 
come more intense inside hotspots when the curvature of 
a NP is more pronounced or when we stack increasingly 
small NPs in what is known as nanolenses. However, 
the volume that enclose these hotspots is reduced at the 
same time and this effect can be so drastic that not even 
a single molecule can fit into the hotspot. With our pro¬ 
posal, we are increasing the electromagnetic fields not 
by sharpening a single NP or using nanolenses, but by 
tuning appropriately the interactions among NPs. An¬ 
other possibility is to use our system for multi-frequency 
simultaneous sensing as it possessed several resonances 
physically separated. 

Mirrored graded-chains can also be used for near field 
spectroscopies without an actual tip. In this case the sys¬ 
tem would act as an effective tip where electromagnetic 
fields are much higher when the analyte is close to the 
center of the array and negligible if the analyte is moved 
away from it. Alternatively it can be used to calibrate 
more standard near-field-spectroscopic devices. 

Its application to solar cells would require a more com¬ 
plex design of the array as one generally need to concen¬ 
trate the electromagnetic fields on a 3D volume. One 
possibility is to use arrays made of multiple 2-D layers 
where the interaction between layers changes gradually 
and the active medium surrounds the NPs that form the 
central layers [16]. Alternatively, one can use the direct 
light-to-heat conversion mechanism mediated by conduc¬ 
tive nano-particles. The dramatic light concentration 
produced by nanostructures such as the one proposed 
here, induces vaporization of the host medium without 
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the requirement of heating the liquid volume to the boil¬ 
ing point [35]. 

Finally, one can also use the proposed structure for 
plasmon-enhanced photo-detectors, helping to reduce the 
size of the photo-detector. This would result in increased 
speed, decreased noise, and reduced power consumption 
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